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In the superclean case the spectrum of vortex core excitations in the presence of disorder is not 
random but consists of two series of equally-spaced levels.B The I-V characteristics of such super- 
conductors displays many interesting phenomena. A series of resonances is predicted at frequencies 
commensurate with the spacing of the vortex excitations. These resonances reveal an even-odd 
anomaly. In the presence of one weak impurity the excitation levels can approach each other and 
almost cross. Absorption at low frequencies is identified with the resonances arising in this case. 
The results of such microscopic theory coincide up to the order of magnitude with both the the- 
ory employing kinetic equationQ and the experiment .□ The non- linear effects associated with Zener 
transitions in such crossings are studied. These phenomena can be used as a probe of vortex core 
excitations. 

PACS numbers: 74.60.Ge,68.10.-m,05.60.+w 



I. INTRODUCTION 

High-T c superconductors (HTSC) in normal state have 
many anomalous properties, differing them from normal 
metals. For example, the relaxation time depends on 
temperature according to the law r _1 ps T/2 (assuming 
h = fes = 1)- It is quite interesting if these properties 
remain anomalous below the transition, when supercon- 
ductivity hinders studying them. Vortex cores retain a 
lot of information about the normal state. However, even 
in the BCS model the properties of vortex cores in super- 
clean superconductors at low temperatures are studied 
insufficiently The purpose of this work is to fill this gap. 

It is believed that the dissipation in the mixed state 
of type II superconductors is associated with the low ly- 
ing excitations arising inside the vortex cores, while they 
are dragged through the sample by the Lorentz forcesB. 
Studying dissipation can therefore shine some light at 
both the spectrum and the relaxation mechanism of these 
excitations. 

Although the excitation spectrum was calculated long 
agcQ, experimentalists apnspached the possibility of its 
observation only recentlyala To this end the excitation 
level spacing luq has to significantly exceed the broaden- 
ing of levels due to various relaxation mechanisms 1/r. 
This superclean case ojqt ^> 1 became possible in Cu 
based high-T c superconductors, as due to a very short 
coherence length £, both ujq ~ l/m£ 2 is large and the 
amount of disorder in the vortex core, contributing to 
the level broadening, is small. It should be noted that 
the measurements in Refs. |3|,|]are done at microwave fre- 
quencies to avoid the effects associated with pinning. 

The DC absorption in superclean superconductorsJi 
been studied theoretically in a number of papers.EJt 
Among them the most notable is Ref. |^, in which the 
excitations were treated in the quasiclassically employ- 
ing the kinetic equation in the r- approximation. The 
result of Ref. for s-wave superconductors is 
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It should be noted however that contemporary high- 
Tc superconductors are extremely anisotropic. For this 
reason the excitations in vortices become quasi two- 
dimensional (2D). It is known that in 2D the kinetic equa- 
tion describes the scattering on impurities poorly. This 
is because a particle has a large probability to return to 
the original position and scatter on the same impurity 
many times. Therefore this problem can be treated i.e. 
considering the transitions between discrete levels. 

Such treatment was realized in Ref. 0. They used the 
random matrix theory to describe the transitions between 
excitation levels. Therefore their theory is not applicable 
to the superclean case. 

The effects associated with the discreteness of the exci- 
tation spectrum in the superclean case have been added 
to consideration in Refs. [l]|| According to these ref- 
erences the discreteness of levels makes the absorption 
non-ohmic. Ref. |6| considers the shake up created by a 
weak impurity passing through the vortex core using the 
Fermi golden rule. They identify a critical velocity below 
which the non-ohmic effects should become pronounced. 
It can be estimated as 



(1.2) 



where k-p is the Fermi momentum. Indeed the impurity 
passing through the core at velocity v can create tran- 
sitions between states separated by — v/Xp, as Ap is 
the characteristic spacial frequency of the wavefunctions 
of the excitations. In the case f2 < lu no transitions can 
occur and the absorption is exponentially small. The au- 
thors of Ref. [I] argue that even if impurity is weak, as 
long as the Born parameter of the impurity 9 3> Af/£, 
due to the special form of the excitation wavefunctions, 
the levels inside the vortex always cross. These crossings 
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occur in the so-called dissipative region, when impurity 
is at the distances a m £8 from the center of the vortex. 
Here the Born parameter < 1. The excitations in this 
case occur due to Landau-Zener transitions happening in 
the level crossings. Due to this effect the absorption at 
v < v c is not e xpo nentially small but is much larger than 
given by Eq. (1.1). 

It should be emphasized that the calculations in 
Refs. have been done for the DC case. The ex- 

periments in Refs. |^,|5| are done at microwave frequen- 
cies (uj k 1AK). The AC absorption in the super- 
clean case has been studied in Refs. |[ [| and [h| All 
these studies are done using the kinetic equation in the 
r-approximation. This paper is dedicated to the micro- 
scopic study of the AC absorption in the superclean case. 
We adopt the mechanism of absorption used in Refs. [l]|], 
i. e. the motion of vortex relative to the impurities brings 
about transitions of the excitations to the higher energy 
levels. 

We study the global diagram of absorption as a func- 
tion of frequency and amplitude of the applied current. 
We find that if the energy relaxation time r £ is large, the 
region on the diagram where kinetic equation gives the 
correct order of magnitude of the result becomes small. 

The outline of our paper is as follows. In Sec. || we re- 
view some basic facts from Ref. [j] about the vortex core 
excitations in the presence of impurities in superclean 
layered superconductors. It is easy to see that in the su- 
perclean case the number of impurities per vortex core 
per crystalline layer can be estimated as A/j ~ I/luqtO 2 . 
If 6 ~ 1 then it is very improbable to have more than one 
impurity inside the core in one layer. Another simplifica- 
tion comes from the fact that if coupling between layers 
is small (open Fermi surface) , the excitation spectrum in 
the presence of impurity can be calculated independently 
for every layer. For this reason Ref. |l| treats the exci- 
tations in the presence of impurity as belonging to one 
two-dimensional layer. As a result they obtain that in the 
presence of impurity the usually equidistant spectrum of 
excitations, pertinent to the two-dimensional clean vor- 
tex core, ceases to be equidistant. However the spectrum 
remains to be strongly correlated. It is shown that the 
system of odd levels and the system of even ones sepa- 
rately continue to be equidistant with the level spacing 
2ojq in each individual subsystem. 



In Sec. IV be describe the resonances occurring in vor- 
tices under the influence of low amplitude, high frequency 
field. The amplitude of vortex motion xq is assumed to 
be much smaller than Af, and frequency of external field 
lu is comparable or larger than ujq- We argue that the 
shape of resonant curves reveals an even-odd anomaly. 
If u> w 2nwo, where n is integer, the transitions occur 
only within each individual subsystem of even or odd 
levels. In this case the resonance is very sharp, with 
the resonant curve determined by the remnant inelastic 
processes. If on the other hand u> ~ (2n + 1) ojq the tran- 
sitions between two subsystems of even and odd levels 
can occur. In this case the resonant frequency depends 



on the position of the impurity, and after averaging over 
this position the resonant curve of absorption becomes 
smeared. 

In Sec. we study the small amplitude low frequency 
absorption. In this case the transition can occur only in 
dissipative regions, where impurity makes even and odd 
levels cross. The result for Rccr^ obtained in this case 



coincides with Eq. (1.1) in the order of magnitude. It is 
therefore purely ohmic. 

The non-linear effects are associated with an increase 
of the amplitude xq. They are of two types. The first 
is attributed to the saturation of energy absorption at 
long timesEi It therefore effectively decreases the mag- 
nitude of energy dissipation. This non-linear effect can 
be neglected if wr e <C 1, where r e is the time of energy 
relaxation. In the latter case another non-linear effect 
becomes important. It arises due to Landau-Zener tran- 
sitions between the crossing even and odd levels, as dis- 
cussed in Ref. [l]. It therefore lead s to an increase of ab- 
sorption with respect to Eq. (1.1). In Sec. VI we present 
a phase diagram of various regimes of dissipation arising 
in this case . 

Sec. VII is dedicated to our conclusions. We discuss 



the possible corrections to our results brought about by 
interlevel coupling, pinning, and d-wave order parameter. 
We compare our results to the existing experiment and 
discuss conditions at which resonances and non-linear ef- 
fects can be observed. 



II. VORTEX CORE EXCITATIONS IN THE 
PRESENCE OF AN IMPURITY. 

In this Section we briefly review some facts about ex- 
citations inside the vortex cose. They can be described 
by the Bogolyubov equations^: 
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Here A (r), /i, V (r), and a are the order parameter, the 
chemical potential, the impurity potential, and the posi- 
tion of the impurity respectively. As it is mentioned in 
the introduction, in the superclean case we can consider 
no more than one impurity per vortex, per layer. We 
will assume that the magnetic field is weak (B <C -ffe) 
and therefore can be neglected in Eq. (2.2). We will also 
assume the s-wave order parameter to have the form 
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where r and <f> are the polar coordinates. 

The low energy excitation spectrum without impurity 
is well known D 



E 



where 
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The corresponding wavefunctions are given by 
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with J ra (x) being the Bessel function and C being the 
normalization constant: 
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If Kramer- Pesch effect takes placeEl at low temperatures 
(T < T c ) A(r) w A(r = oo) = A DO , r>(^. Therefore, 
-ftT(r) is given by 



A 



r 



Consequently 
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The excitation spectrum in the presence of a short 
range impurity at point a has been obtained in Ref. [I]. 
The energy spectrum is given by the following equation 



cos 



nE 



2 ■ tt 2 !/! 2 



4w, 



where 



I = h+ ih = 



2C 2 
irkFa 



,-2K(a) 



V (2k F ) e 



(2.10) 



(2.11) 



and V (q) is the Fourier tr ansfo rm of the impurity poten- 
tial. In the derivation of ( 2.10|) it has been assumed that 
a ^> fcp 1 . 

Note that in the abs ence of impurity (1% = I = 0) the 
spectrum given by Eq. (2.10) is equidistant and coincides 
with Eq. ( |2.4| ). However, in the presence of impurity the 
spectrum ceases to be equidistant. This is illustrated in 
Fig. ^, where the energy levels are shown as functions 
of the distance of the impurity from the center of the 
vortex. 




FIG. 1. The excitation energy levels as functions of the 
distance of the impurity from the vortex center a. The pa- 
rameters used are: Af/£ = 0.1 a nd OAoo/uJokp^ = 0.9 [see 



explanation following Eq. (2.12) below]. Only the vicinity 
of the dissipative region is shown. 

However, the spectrum remains strongly correlated. It 
is easy to see from Eq. (2.1C) and from the Fig. [I] that 
it comprises two series of equidistant levels. The spacing 
within each series is 2wo, while with respect to each other 
they are shifted by a phase depending on the impurity. 

Another feature of the spectrum evident from the Fig- 
ure is that when the impurity is close to the center it 
bring about periodic anticrossings of the levels. The min- 
imum distance between levels in such anticrossings SE, 
according to Ref. ^, has a minimum at point a — ao and 
is determined by the equation: 



SEM^\a-ao\, 
da 

\Kao)\= 2 -^. 

7T 



(2.12) 



Eq. J2.12Q together with Eqs. (gjj) and fl2.1l| ) give 
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Here 9 = mV (2fcp) is the Born parameter. The region of 
the vortex near a = ao is the region where levels approach 
each other very closely. For this reason Zener transitions 
are very probable there. Therefore it was called in Ref. |l| 
the dissipative region. 
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III. MOTION OF THE VORTICES 

Below we recollect a few facts pertinent to the absorp- 
tion by vortices. We consider the system of vortices in an 
alternating electric field oriented in the plane of the lay- 
ers. The magnetic field B is perpendicular to the layers. 
If pinning is negligible the velocity of all of the vortices 
is the same. Let us denote it 



V (t) = v cos (utt) . 



(3.1) 



Therefore the position of impurity relative to the vortex 
is given by 



a = a 



Xq sin (tot) , 
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where Xq is the amplitude of vibrations related to the 
electric field by 
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It is convenient to write the Schroedinger equatio n fo r 
the time-dependent hamiltonian Tt = Tt(t) [see Eq. (|2.2|)] 
in the basis of eigenfunctions of this hamiltonian consid- 
ering time a parameter. These eigenfunctions \n (t)) and 
the corresponding eigenvalues E n (t) can be used to ob- 
tain thet differential equation for the occupation of these 
stateaI3 
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Here 4> n = £E(f)dt'. 

In Sec. [V and [V] we will consider cases when x$ <C Af- 
The transition probability between levels m and n in this 
case can determined by the Fermi golden rule 
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Here the perturbation to the hamiltonian in the system 
of reference moving with the condensate 
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The energy absorption due to one impurity is given by 
Qi (a) = w w m ->„ [f n - f m ] , (3.7) 
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is the Fermi 



distribution function. The absorption averaged over 
many impurities and vortices is determined by 



Q = ni« v J d 2 aQi (a) . 
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Here rii is the three-dimensional (3D) concentration of 
impurities, while n v is the 2D concentration of vortex 
lines. 

Let us now establish the connection between the elec- 
tric field in the plane of the layers £(t) and the velocity 
of the vortices v (t). The AC electric field in a supercon- 
ductor can be written as a sum of two terms: 
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where A is the London penetration depth and j s is the 
supcrcurrcnt. Here the first term is the London electric 
field in the inertial system of reference in which vortices 
are not moving. The other term arises due to the flux 
carried by the vortices with velocity v. 

Let us supplement this formula with the general equa- 
tions for the electric field and dissipation: 



Q = Rep X xj 2 , 



(3.10) 



where the line denotes the time averaging. Further we 
consider the case of small temperatures (T T c ), when 
the relative contribution to the current from normal elec- 
trons is exponentially small. Therefore we can assume 
3 - js- 

In the superclean case Imp xx and p xy are weakly de- 
pendent on the density of impurities. Hence they can 
be fairly well approximated by the values for clean sys- 
tem. Another simplification appears if the spectrum of 
electrons is isotropic^ Then j s — en s v s and A~ 2 = 
47re 2 n s /mc 2 . Eq. ( [3.10 ) can then be rewritten as the 
equation of motion of the center of mass of the conden- 
sate: 
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Here rj is the effective viscosity. Calculation of this viscos- 
ity as a function of frequency and velocity is the purpose 
of this work. Notice that it has both real and imaginary 
components. The real part can be calculated from the 
absorption: 



Q(uj) — Rerjv 2 . 



(3.12) 



The imaginary part is somewhat irrelevant due to the 
fact that it is co mpared to the London term [left-hand 
side of Eq. (3.11)]. At the same time Imry has to vanish in 
the absence of impurities as in the opposite case it would 
shift the position of_the cyclotron resonance, violating 
the Kohn's theorenO. Therefore in the superclea n case 
it is small compared to the left-hand side of Eq. (3.11) 
and can be neglected. 



Eq. ([T9|) together with Eq. ( |3.11| ) result in the follow- 
ing equation for the vortex motion: 
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Here lo c = \e\B/mc is the cyclotron frequency. In the 
superclean case the ratio r)/muj c is proportional to the 
density of impurities. In the majority of cases considered 
below this quantity is small. For example the kinetic 
equation at zero frequency gives the following estimate 
for this quantity: 



ri/muj c — 1/loqt <C 1. 



(3.14) 



Our calculations discussed below confirm this result. 
Therefore for the superclean c ase it is plausible to assume 
that v is close to v s [see Eq. ( 3.13 )], i.e. the condensate 
and the vortices almost do not move with respect to each 
other. 



Eq. ( 3.11 ) can also be used to obtain an expression for 
the resistivity tensor: 
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The dissipative component of the conductivity tensor fol- 
lowing from this expression is given by: 
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where 7 = Herj/mtOc is smal l in the superclean case 
in a ccord ance with Eq. ( |3.14| ). At small frequencies 
Eq. (3.16) can be seen to describe the cyclotron reso- 
nance. The cyclotron resonance in superconductors was 
observed experimentally in Ref. [l5] and was described 
by Ref. [lO] in the framework of the kinetic equation in 
the r-approximation. Note that the cyclotron resonance 
rendered by this expression is sharp if 7 <C 1 or ex- 
actly in the con ditions of superclean case. Note also 
that as Eq. ( 3.16 ) is a conseque nce o f quite general equa- 
tion of moti on of condensate (3.11). In the derivation 
of Eq. (3.16) we have neglected by the imaginary part 
of r\ in comparison to rruo. It is possible to do so be- 
cause Imry has to vanish in the absence of impurities. In 
the opposite case it would shift the position ofi-the cy- 
clotron resonance, violating the Kohn's theoremliJ, as it 
was mentioned before. 

Let us analize the influence of pinning on the motion 
of the condensate. We will ignore the dissipation for a 
moment. If pinning is present the pinning force has to 
be added to the total force acting on the condensate in 



Eq. (3.11) 



F p = —iav/uj. 
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Here iv/ui is the displacement of the vortex lattice (as- 
suming it to be rigid), a is the pinning parameter, related 
to the critical current j c by 
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In the last expression we have introduced the pinning fre- 
quency Up. This frequency can be related to the critical 
velocity v p = j c /n- s \e\ and the cyclotron frequency by 



(3.19) 



Eqs. (3J3) and ( 3.11 ) give the following expression for the 
modified by pinning frequency of the cyclotron resonance: 

Co = — £ (3.20) 



At a non-zero ui p the 



In the limit co p — > 0, Co — > uo 
cyclotron resonance occurs at a frequency larger then uo, 
(Ref. [h]). At t he same time the equation of vortex mo- 
tion, Eq. (§t|) 
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Therefore the pinning is important if uo <C uo 2 /uo c . 

Pinning can be ignored at large frequencies (uj ^ > Co) or 
if velocity is larger than j c /N s \e\. Then Eq. ( 3.16 ) gives 
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Using this expression one can relate the resonance struc- 
ture arising in Re/? xx near to ~ loq to that in RetTxx. This 
resonance is described in the next Section. 



IV. RESONANCES AT HIGH FREQUENCY AND 
SMALL AMPLITUDE. 

In this Section we will assume the displacement of the 
vortex due to the driving field to be much smaller than 
the smallest scale at hand Af- In this case to calculate 
absorption it is only natural to employ the Fermi golden 
rule. The absorption averaged over a large interval of 
lo in this linear response formalism is of the same or- 
der of magnitude as given by Eq. (1.1). However at the 



frequencies equal to multiples of loq resonances occur in 
the sample. The aim of this Section is to calculate the 
absorpt ion near these resonances. 

Eq. (|3.5| ) emphasizes that at the conditions of reso- 
nance the energies of two states E n and E m should be 
close to the multiples of loq. This implies that they should 
not be disturbed strongly fr om t he values in the absence 
of impurity given by Eq. ( |2.4| ). Therefore near reso- 
nances the main contribution to the integral in Eq. ( |3.8| ) 
comes from the regions far from the center of the vortex 

~ 0-0 <C a. In these regions one can treat the influence 
of impurity perturbatively. In the first order of the im- 
purity potential we obtain the following expressio ns fo r 
the correction to the energy of excitation from Eq. ( 2.10| ) 
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(4.1) 



At these conditions the transition probability in the first 
non-vanishing order in the impurity potential can be de- 
termined taking as wavefunctions the sta tes o f the hamil- 
tonial with no impurity [Eq. (Bjj)]. Eq. (3.6) then gives 
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Here <f> is t he a ngle between a and xq. The golden rule 
expression (3.5) can be readily rewritten as follows 
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where Aw = w — Iloq < wo/2, with I being integer, is the 
deviation of the frequency from the resonance, and index 
of the m atrix element should be chosen in accordance 
with Eq. 

Looking at the <5-function in this expression one imme- 
diately sees that the answer should be quite different for 
w close to even and odd multiples of Wo (for even and odd 
I). Indeed if I is even the resonant 5- function reduces to 
S(Au). Therefore the condition of resonance is the same 
at any position of impurity relative to the center of vor- 
tex. Hence 6 (Aw) appears as a factor in front of the 
expression for absorption. Therefore the resonances are 
very sharp in this case. 

In the opposite case, when I is odd, the resonances can 
occur even if Aw is not zero, due to the shift of energy 
levels by the impurity. Thus for odd I the resonances are 
broadened by the presence of impurity. This difference 
between even and odd I is a consequence of the men- 
tioned in Sec. || property of the spectrum. The series 
of odd and even levels are shifted with respect to each 
other by the impurity, preserving equal spacing of 2wo 
within each series. Thus the resonances at even frequen- 
cies occur within each series of levels and are sharp, while 
the perturbation at odd frequencies mixes even and odd 
levels, making resonance broader. Below we study these 
two cases separately. 

A ) Odd to /luq . 

Eq. (|4.3| ) for this case can be rewritten as follows 
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The resonances in the integrand of the last expression 
occur very frequently (once in Af/4). One can average 
contribution from these resonances over much larger in- 
tervals to obtain an integral over a slowly varying inte- 
grand: 
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The integral can be evaluated by changing the vari able o f 
integration from d to \I\. They are related by Eq. ( 2.11 ). 
The function a (|/|) can be found asymptotically in two 
limiting cases 



a(\I\)={ 



2C 2 
Trk-p\I\m ' 



\I\ > 0u>o 



Here 



mV (2k F ) < 1 
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is the Born parameter of the impurity. Calculating the 
absorption in these limiting cases one obtains 
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for Aw (9wo. It is convenient to express this answer in 
the units of absorption obtained from the kinetic equa- 
tion. Assume that r is determined by scattering on short 
range impurities: 
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with v = raj'K bein g 2D density of states in the layers 
Then from Eq. ([0]) one obtains 

Qk = tt^TT— = «i^ 2 -. (4.11) 

2n s |e|cwoT wo 

To simplify our consideration we wi ll ac cept t he K ramer- 
Pesch ansatz, expressed by Eqs. (2.8) and (|2.9|), after 
what we obtain 
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for 9u)q <S Aw <C wo and 
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for Aw <C Oujq. 

It can be noticed from these equations that when the 
frequency approaches the resonance the absorption first 
increases in accordance with Eg. ( 4.12 ) and then de- 
creases according to Eq. (4.12). The former expression 
therefore describes resonance near wq, while the latter 
describes the antiresonance. The numerical evaluation 



of Eq. (4.5) shows that the absorption reaches its maxi- 
mum of max(0.2/6»; 0.6) at point Auj w 0.59w . K > 0.3 
the resonant part is not very well pronounced. On the 
other hand the antiresonace at Aw < 0.59wq exists in- 
dependently of the value of the Born parameter. An an- 
tiresonant behavior has been described earlier in Ref. [l^. 
The width of the antiresonance obtained in this reference 
is 1/t. We obtain for the width of the antiresonance 
Ao; « 0.59loq. The discrepancy is the consequence of the 
mentioned failure of the r-approximation in (quasi) two 
dimensions. 

Another perturbation causing the resonance could be 
due to the inertial forces. This perturbation is impor- 
tant because it exists even in the cores with no impuri- 
ties. The co rres ponding correction to the hamiltonian, 
similarly to (|3.6|), is 



5Hm 



pv s ; a; VA 
aj VA*; pv s 



(4-14) 



Here the diagonal terms represent the Doppler shift of 
the energy of the quaziparticles caused by the condensate 
moving with velocity v s (Ref. while the off-diagonal 
terms are associated with the motion of vortex itself.u 



The matrix element of ( 4.14 ) calculated between the pure 
states relevant to the neighborhood of the resonance can 
be rewritten as follows: 



(m\ SHincrt \n) = (to | p (v s - v) cr \n) 



(4.15) 



Here <jq is the unit Pauli matrix. In this equation we 
used the identity: 

(m\ cti^oVA \n) = (m| [V; H] \n) = (m\ pva \n) . 

(4.16) 

The absorption due to this perturbation is given by: 



Q 



Inert 



^0 (V - VsY 

' 4r £ Aoj 2 + Ts 



(4.17) 



where v — v s is given by either Eq. (3. IS) or by Eq. (3.21). 
In case if pinning is ne gligible v s is very close to v, as it 

from 
mertial 

forces can then be estimated as 
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follo ws fr om Eq. ( |3.13|) . Taking an estimate for r\ 
Eq. ( 4.13| ) at Auj ~ r" 1 the absorption due to the ins 



Qlnert/QK ~ l/^o rT e- 



(4.18) 



This is just a small correction to Eq. (4.13). The signifi- 
cant reduction of the absorption is a consequence of the 
Galellian invariance. 



B) Even uj/loq. 

In this case due to ho moge neous broadening of the lev- 
els the 8 function in Eq. (4.3) is replaced by the following 
expression: 



5(Auj) 



1 



7T7> ( Auj 2 



(4.19) 



Eq. ( 2.11 ) gives the following result for the absorption 



4C 4 6> 2 
kn v — ^ — 



■8 (Aw) In 



(4.20) 



Using the Kramer-Pesch ansatz (2.8) and (2.9) one read- 
ily obtains 



In 



Q 



Qk ln ^ 

T 



8(Aw) 



(4.21) 



Note, that the absorption averaged over frequency for 
w > loq is the same in both even and odd cases and 
equals to 



1 M+l/2)u l n 
<?K <2kM) J(l-l/2)uo ^ ln^ — 

T 



(4.22) 



V. SMALL AMPLITUDE LOW FREQUENCY 
ABSORPTION. 

At a frequency much smaller than loq absorption is only 
possible if the impurity makes two levels approach each 
other to a distance smaller than w. This is possible if 
impurity is situated within dissipative region described 
in Sec. Q). From Eq. ( |2.12| ) it follows that the condi- 
tion SE < w, lu — > can be satisfied in a narrow region 
adjacent to a = ao- In this region odd and even levels 
periodically anticross, approaching each other to small 
distances. 

When two levels are anomalously close all others can 
be ignored and the two-state system can be described by 
a hamiltonian 



H = 



8 —a 



(5.1) 



where 8 > is a slowly varying function of the a, while a 
changes rapidly and a — at the point of level anticross- 
ing. The eigenvalues of this two-level system are given 
by 



(5.2) 
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Comparing this equation with Eq. (2.10) in the neigh- 
borhood of the anticrossing one obtains the following ex- 
pression for a and S 



S = 



SE 
h p 



2kpSaujQ 



(5.3) 



where SE is given by Eq. (2.12) and 5a is the distance 
from the point of anticrossing. The cigenmodes of the 
hamiltonian are described by 



/ l\e\±6 \ 



a± 



2\e\ 



± 



\s\tS 



(5.4) 



V V 2 N / 



The upper and lower signs in this expression perta in t o 
the states with positive and negative energy in Eq. (5.2). 
Now we c an apply the golden rule expression for absorp- 
tion (|3.8|) to this two-level system: 



(5.5) 



Q = nin v oj J d(f)ada—x Q cos 



dH 
da 



5[lo - 2s (a)], 



with 



— -) = -k ¥ u — 

oa 7T £ 



(5.6) 



One can average over frequently occurring resonances in 
the integrand of Eq. ( |5.5| ) to obtain an integral over a 
slowly varying function 



Q = 2nin v XQk ¥ uJoaQ 
S 2 (a) 



da 



(5.7) 



Ch angin g the variables of integration from a to S using 
Eq. 2.12, i.e. using 



da — AdSda/d\I\ = 2naodS/uJo, 

we eventually obtain 

Q _ 
Qk 



K, In ( ^ 

T 



(5.8) 



(5.9) 



To obtain thi s an swer we have employed the Kramer- 
Pesch ansatz (O) and (Oh. 



Eq. ( |3.8| ) can also be evaluated for general u. The 
result reads 



Q 



2uj 



2fioor 



2^r 



LJ 



Sill 



(5.10) 



where floor(cc) is the maximum integer smaller or equal 
than x. Note that both asymptotics given by Eqs. ( 4.12| ) 
and (5.G) can be obtained from this expression in the lim- 
its lo — > and lj — > (2n + l)yp, n = 0, 1, . . . respectively. 
This result is shown in Fig. 




FIG. 2. Abs orpti on as a function of frequency deter- 
mined by Eq. (5.1C). At add frequencies only resonant 
behavior is sho wn. The sharp resonances at even frequen- 
cies [Eq. (4.21)] are shown by vertical lines. The non-zero 
limitQ((jj — > 0) is evaluated in this Section. 

As it was mentioned in the previous Section, due 
to numerical reasons, the absorption at odd frequen- 
cies (uio, 3u>o, . . .) does not raise above max(O.2/0; 0.6). 
The resonances at odd frequencies can therefore be pro- 
nounced on ly if 9 < 0.3. We expect that in practice 
9 < 1. Eq. ( 5.10|) and Fig. ||can therefore be understood 
only asymptotically in the limit 9<1. 



VI. NON-LINEAR EFFECTS IN THE 
ABSORPTION 



The non-linear effects in the absorption are interesting 
since in many cases they depend on the energy relax- 
ation time of excitations. Therefore they can be used for 
studying this quantity. We attribute the non-linearities 
to two different phenomena. The first phenomenon is 
the saturation of absorption due to the redistribution of 
level occupations. It corresponds to the deviations from 
the golden rule expression ( |3.7| ) . This leads to a decrease 
of the dissipative component of resistivity with respect 
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to a linear response result. The second group of non- 
linear phenomena includes the effects associated with di- 
rect transitions between levels due to their anticrossings. 
They increase the absorption. 

Consider the deviations from the golden rule expres- 
sion first. In a two- level system the resonance with ex- 
ternal field brings about the rotation of the population 
of the levels with Rabi frequency Or. For the cases de- 
scribed in Sec. IV and [v| the Rabi frequency is given by: 



fl R = \(m\SH\n)\ 
x 



(6.1) 



LU 



At 



If the perturbation is applied during a time interval 
longer than 1/f^R, the absorption saturates. In this case 
the system can absorb only if there are inelastic relax- 
ation mechanisms present. If the corresponding relax- 
ation time is r E , the energy u>, corresponding to the ex- 
ternal field frequency, is absorbed once in r e . This is in 
contrast to the golden rule expression, from which this 



time is seen to be I/Or. To account for this this effect 
the absorption should be renormalized as follows 



Q 



Q 



(6.2) 



The denominator in this expression is large if v/v c 3> 
uj/ujqT £ . The absorption in this case is 



Q~Qk 



v c uj 



(6.3) 



The non- linear effects are relatively abundant. To visu- 
alize them we draw a diagram in the space of parameters 
v and u>. This diagram is shown in Fig. |^. It displays 
the absor ption measured in the units of Qk given by 
Eq. (pd| ). 

The golden rule absorption considered in the previous 
two Sections is situated in sectors (a) and (b) of Fig. g. 
The saturation of the golden rule express ion becomes im- 
portant on the line v/v c = uj/uJ^T e . Eq. ( |6 - 3| ) determines 
the absorption in sector (d). 



v/v c 



1/(co„t e ) 2 




FIG. 3. The diagram of absorption. The values of 
Q/Qk are mapped for different regions in (v/v c , ui) plane. 
The equations for the boundaries between these region are 
made italic for convenience. 



We now turn to the discussion of the influence of 
the Zener transitions between anticrossing levels. These 
transitions occur within the dissipative region, where 
some levels can be anomalously close. In this case the 



system of two almost degene rate states can be described 
by a two- level hamiltonian (5.1), and their energies are 
given by Eq. (5.2). The occupations of such a two-level 
system before the anticrossing (pi,j?2) and after the an- 
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ticrossing (p^Pg) are related as follows 



P2 



1-W; W \(pi 
W- 1 - W )\ P 2 



(6.4) 



The transition probability to the other level relevant for 
the absorption is 



W = exp 



tt 2 SE 2 



(6.5) 



Here SE is th e sm allest distance between levels deter- 
mined by Eq. ( [2.12 ). 

Zener transitions can occur if W ~ 1. This determines 
the width of the dissipative region as well as the number 
of effectively working anticrossings there N. Indeed, if 
an impurity is situated in a ring of width XpN and ra- 
dius ao ~ 9£ it can produce Zener transitions. We obtain 
from Eq. flfl^ ) 



SE 



4d 

da 



Then from Eq. (ETq) we have 



(6.6) 



(6.7) 



One can distinguish two cases, depending on whether 
N is large or small. If N is large, an impurity passing 
through the dissipative region brings about many tran- 
sitions. The opposite case N <C 1 can be realized if 
velocity is small v <C v c (Af /0£) 2 . Below we first consider 
this latter case in some detail. 

If N <C 1, for a typical impurity Zener transitions are 
improbable. However if the Born parameter of the impu- 
rities has a little dispersion S9 = J (9 2 ) — (9) 2 ^ 0, the 
Zener transitions still can happen, yet on a very small 
fraction of impurities. Indeed, the current position of 
the level crossing a is quantized in the units of Ap/4. 
The uncertainty of ao ~ Ot; is Sao ~ 59^. Let us look at 
Eq. ( 2.12j ). If So,q > Ap/4, SE in this expression can acci- 
dentally become zero. This condition can also be rewrit- 
ten as 89 > Af/£ <C 1. This assumption can therefore be 
realized even if the uncertainty of the Born parameter is 
small. 

This can be reformulated in terms of an assumption 
about the distribution function of SE. In the case when 
56 > Af/£ this distribution function p(SE) has a non- 
zero limit when SE goes to zero. In particular, taking 
Eq. (|6.6[) into account, we obtain 



p(6E) 



dN 



d[8E\ Apc^o 



The expression for the total absorption Eq. (|3.8| ) can be 
rewritten in this case as follows (a ~ 6*£, da ~ xq) 

/>OC 

Q~ ni n v 6€x Q 1 (6E)p(SE)d[SE], (6.9) 
Jo 



where Q\(8E) is the absorption due to one impurity driv- 
ing Zener transition on a couple of levels with the shortest 
separation SE. Our task from now on will be to estimate 

Qi- 

For the following discussion it is convenient to consider 
two cases xq <C Ap and xq ^> Ap separately. In the first 
case an impurity can cause only one Zener transition per 
period, while in the second case there can be many. 

A) Xq < Ap. 

Consider an impurity at a distance less than xq from 
one of the anticrossings. Such an impurity can cause 
Zener transitions once in uj~ 1 . The absorbed energy, 
however, depends on the relation between the charac- 
teristic time of Zener transition At and the inelastic re- 
laxation time 7>. The former can be estimated as 



SEX F 
At - 

UJQ v 



(6.10) 



and is the time during which two levels are at a distance 
of the order of SE. The absorption rate due to one such 
impurity can be estimated in three limiting cases as fol- 
lows 



f 5EWu;, t e < At 



, , I -r^UJoWuj, At < T e < oj- 1 

yi ~ < Ap 



2^0 1 i 

— lu W— , Ai<w 1 <C r e 

Af T e 



(6.11) 



As it has been mentioned above this is true if the impu- 
rity is at a distance smaller than xo from the anticrossing. 
In the opposite case the Zener transitions cannot occur 
and the absorption is negligible. Eqs. ( |6.9| ) and (6.11) 
then result in three limiting cases 



UJ T e 



V 1 
UJT e < 1, — > 



Q - Qk < 



,! 2 T 2 

J T e 



V 1 
UJT t < 1, — < -5-2 (6.12) 

Vc T e 



U! 2 U!oT e V V, 



U)T e » 1. 



These limiting cases correspond to the regions labeled (f ) , 
(c), and (e) (Af xq) correspondingly in the diagram 
in Fig. § 

Note that for this case the answer does not change 



(6.8) when velocity passes across the value 



where N ~ 1. This is due to the fact that if xq <C Af even 
for the case N S> 1 the moving impurity cannot produce 
more than one Zener transition per period 2it/uj. This 
happens not to be the case for xq S> Af , when the bound- 
ary v = u c A F /9 2 ^ 2 does exist. 
B) x > A F . 
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As it has been mentioned above for this case the divi- 
sion into small and large velocities is essential. We start 
from small velocities v <C v c Xp /9 2 £ 2 , as in the experi- 
mental conditions^ this is the region where the non-ohmic 
absorption can first be observed. 

For small velocities a consideration analogous to 



Eq. (6.11) leads to 



SEWlo, 



7v < At 



Qi ~ tanh(wr £ )< Af 

LUT £ 

^-luqW— , A* < Af/u < r £ 
V. A F T e 

(6.13) 

The factor tanh (wt e ) /wt e accounts for the saturation 
of absorption when u> 3> t~ . This is analogous to the 
saturation discussed in the Fermi golden rule case. This 
results in 



Q ~ — ^ tanh (wr £ ) < 

W7V 



— , 1/W T £ < U 



« 1/22 , / (6-14) 



1. 



« « l/^r e 2 



Note that for small frequencies lj -C l/r e we have repro- 
duced two of the results from the previous case io < Af, 
namely, sectors (f) and (c) of the diagram in Fig. [|. 
For this reason the line xq = Af in this Figure is made 
dashed, as it does not divide physically different regions. 
We have also obtained a new result for sector (g). 

Let us now discuss the case N 3> 1. If the amplitude 
of vortex motion xq exceeds XpN, the impurity can cause 
N Zener transitions per period of motion 2tt/uj. This im- 
plies that the excitations can ballistically propagate up in 
energy to the hight ujqN above the Fermi level. There- 
fore the power absorbed in one vortex can be written as 
follows 



Qi ~ — tanh (ut e ) ujqN . 

TV 



Hence, for the case XpN <C x Eq. 
(a ~ 6£, da ~ Xq) 



(6.15) 
|3.8|) results in 



Qk 



A WT £ < 1 

Af 



(6.16) 



■, WT e > 1. 



These answers correspond to the sectors (h) and (i) in 
Fig. | 



The equality xq = ApiV takes place on the line v/v c — 
(ui9£/uioAf) 2 . This line is the boundary between sectors 
(i) and (1) of the diagram. In sector (1) xq <C Af./V, and 
the maximum energy above the Fermi level that excita- 
tions can reach is wq^o/Af- Therefore 



Qi ~ — tanh (lut £ ) uj ( t^- 

1~s \Af 



(6.17) 



Eq. ( |3.8D results in (a ~ 0£, da - A F A) 
Q ~ niriv^ApiVQi 



Qk 



UJT £ < 1 



(6.18) 



V 7 v/v c /ujT e , U!T £ > 1. 



These answers correspond to the sector (1) of Fig. ||. 

Finally we would like to discuss the crossover between 
the regions (h) , (i) , and (1) and the non-ohmic regimes at 
high velocities in the sectors (m) and (n) . When velocity 
exceeds the critical velocity v c the number of anticross- 
ing is restricted by the size of the dissipative region. It 
is therefore 



N 



0L 

A F 



(6.19) 



Therefore the answers in sectors (m) and (n) can be ob- 
tained by renormalization those in sectors (h) and (i) by 
the factor v c /v. 

In the case of small frequencies lut £ <C 1 one can give 
the following interpolation formula for the absorption 



Q_ 

Qk 



1 



^/v c /v/u> t s + \Jvjv c 

1 1 



(6.20) 



yJv c /v/uj Q T £ + 6>£/Ap 1 + v/v c 



The second term arises due to the Landau-Zener tran- 
sitions on rare impurities, causing very small 6E. Such 
impurities can therefore cause no more than one transi- 
tion per period. The third term is associated with the 
cascade of Landau-Zener transitions. If Uqt £ < #£/Ap, 
then the third term is larger than the second and sectors 
(g) and (k) in Fig. || do not exist. 



VII. SUMMARY AND CONCLUSIONS 

In this work we have studied the influence of the dis- 
creteness of excitation levels on microwave absorption in 
superclean layered superconductors. At low amplitudes 
and low frequencies (w <C ujq) the absorption coincides 
with the result of Ref. ||. With increasing amplitude of 
vortex motion two non-linear effects are observed. In the 
case lot £ > 1 an increase of the amplitude decreases the 
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dissipative component of resistivity. This is due to the 
deviation of the occupation of the excitation levels from 
the Fermi distribution. In the opposite case ujt £ <C 1 an 
increase of the amplitude brings about an increase of the 
dissipative component of resistivity due to Zener transi- 
tions. Therefore this effect can be used to measure the 
inelastic relaxation time. 
In Sec. 



Ill 



we discuss the conditions at which the 
cyclotron resonance can be observed in superconduc- 
tors. We find that the broadening of the resonant curve 
is uj c /ojqt. We conclude therefore that the cyclotron 
resonance is sharp if the condition of superclean case 
(luqt 1) is satisfied. This provides and independent 
way to measure this quantity. We also calculate the cor- 
rection to the cyclotron frequency brought about by pin- 
ning. 

For small amplitudes and large frequencies we obtain 
the series of antiresonances at the odd frequencies com- 
mensurate with the excitation level spacing [(2n + 1)uj q , 
n = 0, 1, 2, . . .], and the series of resonances at even 
frequencies commensurate with the level spacing [2nw ]. 
Therefore the resonant behavior of the vortex cores in su- 
perclean superconductors reveals an even - odd anomaly. 
The existence of antiresonances at odd frequencies is as- 
sociated with the retraction of the effectively working 
near the resonance impurity from the vortex to the dis- 
tances a>^. There it creates very small matrix element 
of transition between vortex states. On the other hand 
at even frequencies the effectively working near the reso- 
nance impurity resides inside the vortex core (a < £). 

Let us discuss the assumptions that we have made. We 
assumed the temperature to be sufficiently small. First, 
it has to be much smaller than T c , so that the relevant 
excitations are given by Caroli - de Gennes - Matrison 
theorvQ. Second, r e (T) should be large enough to satisfy 
the condition of the superclean case ujor e 3> 1. 

We also adopted the model of disorder consisting of 
strong short-range impurities with the Born parameter 
satisfying the condition 8 3> (uiot)~ 1/>2 . At this condi- 
tion there is no more than one impurity per vortex core 
per layer. In the opposite case of white noise disorder 
potential t he ex citation levels can be shown to be broad- 
ened by yjojo/r. This is analogous to the case of Landau 
level broadening by the white noise disorder considered in 
Rcf. [l7| In this case we expect the absorption due to im- 
purities. At low frequencies (to Cwo) to be exponentially 
small.taij 

We have considered 2D case assuming that tunneling 
between layers is small. In the presence of such tunnel- 
ing the excitation levels with no impurities are broad- 
ened into band. The width of the band is of the order 
of 5wo ~ ujoe, where e -C 1 is the anisotropy parame- 
ter. However, in the presence of impurity one level can 
leave the band and become discrete. The virtual tran- 
sition to the other layers produce only small correction 
to the energy of this level and to our results. This cor- 
rection can be accounted for by the renormalization of 
the impurity potential. However the antiresonances at 



(2n + l)<x>o and resonances at 2nuo should be smeared 
in the presence of anisotropy. The finite homogeneous 
broadening of levels r~ 1 should also add some broadening 
to the (anti)resonances. We therefore conclude that the 
(anti)resonances should be broadened by min(<5o;o, T e -1 )- 
We disregarded pinning of the vortex lattice by the 
impurities in our consideration. It is possible to do so if 
the pinning does not affect motion of the lattice. This 
happens when ui 2 < ujlo c , where uj p is t he c harac terist ic 



pinning frequency introduced in Section [Hl| [Eq. (3.20)] 
Pinning can also be neglected if the current significantly 
exceeds the critical current. 

We have assumed the s-wave pairing mechanism in our 
treatment. In practice the superclean case can be realized 
in high-T c and organic superconductors. It is accepted to 
think that these materials have a d-wave coupling mech- 
anism. Therefore it is important to study the absorption 
in d-wave superclean superconductors. This study has 
been done by Kopnin and VoloviktJ in the framework of 
kinetic equation in the r-approximation. We expect how- 
ever that the differences between the results of kinetic 
equation and the microscopic approach studied in this 
work persist in the d-wave superconductors. Therefore 
we assume that some qualitative results obtained in our 
paper are applicable to these materials also. This matter 
requires a further study. 

One of the known to us experiments on superclean 
samples is described in Refs. [3[ This paper reports a 
large ojqt « 14 observed in 90K single-crystal YBCO 
sample. This quantity weakly depends on temperature 
below 17K. Therefore we expect that at the conditions 
of weak temperature dependence scattering on impurities 
provides the main mechanism of absorption. 

As it follows from our results the study of the frequency 
dependence and/or non- linear effects of absorption it is 
possible to determine the inelastic scattering time r E even 
if it is much larger then the elastic one. Knowledge of 
r e is important for undesstpading of the peculiarities of 
high-r c superconductors 
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